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By using the composite many-body theory for Frenkel excitons we have recently developed, we 
here derive the ground state energy of A'^ Frenkel excitons in the Born approximation through the 
Hamiltonian mean value in a state made of A*' identical Q = excitons. While this quantity reads 
as a density expansion in the case of Wannier excitons, due to many-body effects induced by fermion 
exchanges between A*' composite particles, we show that the Hamiltonian mean value for Frenkel 
excitons only contains a first order term in density, just as for elementary bosons. Such a simple 
result comes from a subtle balance, difficult to guess a priori, between fermion exchanges for two 
or more Frenkel excitons appearing in Coulomb term and the ones appearing in the A'' exciton 
normalization factor - the cancellation being exact within terms in 1/A^s where Ns is the number of 
atomic sites in the sample. This result could make us naively believe that, due to the tight binding 
approximation on which Frenkel excitons are based, these excitons are just bare elementary bosons 
while their composite nature definitely appears at various stages in the precise calculation of the 
Hamiltonian mean value. 
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I. INTRODUCTION 

Quite recently, we have constructed a many-body theory for Frenkel excitons^'^ similar to the one we have developed 
for Wannier excitons. In these theories, the composite nature of the particles is treated exactly through a set of "Pauli 
scatterings" for carrier exchanges between two excitons in the absence of carrier interaction. Exchanges between two or 
more excitons, which control the many-body physics of these composite particles, are visualized and readily calculated 
by using a new diagrammatic representation, called "Shiva diagrams" in relation with their multiarm structure. 
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Alternative approaches which allow one to take into account the Pauli exclusion principle in the case of Frenkel 
excitons are based on expressing the fermionic Hamiltonian as an infinite series expansion in creation and destruction 
operators for excitations on atomic sites. In the formalism proposed by Agranovich and coworkers^ili^, these operators 
obey bosonic commutation rules, while in the method developed by Mukamel and coworkers^i^i^ i ' ^ ^ , the commutation 
rules for these operators are not bosonic. The precise forms of the resulting Hamiltonian expansions produced by 
these two approaches also differ from each other. In contrast to these approaches, the composite boson many-body 
theory we have constructed^^ avoids the difficulty of deriving Hamiltonian infinite series expansions with an increasing 
number of operatorsi^-, the processes involving any number of excitons being in our formalism treated on the same 
footing. This enables us, for the first time in our best knowledge, to self-consistently study interactions between an 
arbitrary number of Frenkel excitons while including the tricky consequence of the Pauli exclusion principle in an 
exact way. 

The composite nature of the Wannier excitons turns out to be crucial in their many-body physics. This becomes 
quite reasonable once we note that the extension of the electron-hole relative motion in an exciton - which is expected 
to control carrier exchanges - is as large as the scale for Coulomb interaction between excitons, this scale being the 
exciton Bohr radius calculated with the semiconductor dielectric constant and the conduction and valence band 
effective masses. In contrast, since Frenkel excitons are made of electrons and holes on the same atomic site, we 
could think that, as the extension of their electron-hole relative motion is essentially zero, fermion exchanges between 
Frenkel excitons should reduce to zero; so that these excitons should behave as elementary bosons. In our previous 
work on the Frenkel exciton many-body theory^, we have shown that, indeed, the closure relation for N Frenkel 
excitons has the same prefactor as the one of N elementary bosons, while for Wannier excitons, this prefactor 

is which is the signature of the fact that Wannier excitons are made of two independent fermions. However, 

for many other basic properties, Frenkel excitons are composite objects definitely. In particular, they do have Pauli 
scatterings for fermion exchange in the absence of fermion interactions. These scatterings lead to a normalization 
factor for N identical excitons different from its A^! value for A^ elementary bosons, by a factor Fn which decreases 
exponentially with A^. This normalization factor thus exibits the same "moth eaten" effect induced by Pauli exclusion 
as the one found for Wannier excitons. 

Fermion exchanges between N Frenkel excitons are found to be controlled by the dimensionless parameter 

V = N/Ns (1.1) 
where Ns is the number of atomic sites in the sample. It is of interest to note that this parameter is actually similar to 



the dimensionless parameter rj = N [a^/ L)^ which controls fermion exchanges between N Wannier excitons. Indeed, 
(L/ax)^ , hke N^, is the maximum number of excitons a sample of volume can accommodate. Let us stress that, 
while the expression of rj found for Wannier excitons could make us think that the parameter which controls fermion 
exchanges is linked to the extension of the electron-hole relative motion wave function, the expression of rj for 
Frenkel excitons rules out this physical understanding - otherwise rj for these excitons would reduce to zero due to 
the tight- binding approximation underlying the Frenkel exciton concept. 

In the present work on Frenkel excitons, we are interested in the density dependence of the N ground state exciton 
energy in the Born approximation. We reach this quantity through a procedure similar to the one we have used for 
Wannier excitons, namely we calculate the Hamiltonian mean value in a N ground state exciton state with momentum 
Q = 

_(MBlHBlW)_ .... 
{v\B^B'q \v) 

where creates a Frenkel exciton with momentum Q. Let us recall that, while the interacting part of composite 
quantum particles cannot be written as a potential so that a "zero order state in the exciton interaction" cannot be 
properly defined, the state B^^ \v) plays this role since {H — NEo)B^'^ \v) cancels if we drop all Coulomb scatterings 
between excitons. As H contains Coulomb interaction at first order by construction, {H)j^ thus corresponds to the 
ground state energy of N Frenkel excitons at first order in Coulomb interaction, i.e., in the Born approximation. 
We are going to show that, while the N exciton normalization factor reads as 

{v\ B^Bl"" \v) = NIFn (L3) 

where Fn, calculated in our general paper on the many-body theory of Frenkel excitonsS', has an expansion in 77 which 
comes from fermion exchanges and makes it exponentially decrease with N, the Hamiltonian mean value is found to 
be exactly equal to 

(^).- ^^0 + ^^,0 (1.4) 

7o comes from the Coulomb scattering of just two excitons due to both direct Coulomb processes and indirect Coulomb 
processes, i.e., the so-called " electron- hole exchange" which is the only process responsible for the excitation transfer 
in the case of Frenkel excitons. Eq. (1.4) thus leads to 

{H)^^N{Eo + riAo) (L5) 



Ao = 70 [1/2 + 0{1/N) + 0{1/N,)] 
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This shows that the Hamiltonian mean value depends on the dimensionless parameter r] associated to exciton density 
through a first order term only, corrections being in N^^ and N^^, but not in N/Ns = rj. 

This has to be contrasted with the Hamiltonian mean value for N Wannier excitons that we have found to expand 



in 3D systems for which the exciton ground state energy is i?o = — -Ro where i?o = /ie^/2£^ ~ 1/2/ia^ is the exciton 
Rydberg. We have explained the presence of terms in density higher than by saying that, while (H) contains 
one 2x2 Coulomb process only by construction, interactions between more than two excitons exist through fermion 
exchanges induced by the composite nature of the particles. 

If we now return to Frenkel excitons, the reason for the disappearance of high order terms is far from trivial. While 
the tight binding approximation underlying the Frenkel exciton concept makes most Coulomb exchange processes 
between excitons reducing to zero, this is not enough to explain the result since the denominator of (H) contains 
a density expansion through Fm- The fact that (H) does not contain terms in 77" with n >2 actually means that 
similar rj" terms exist in the Coulomb part of (H) pj in order to cancel the 77" terms coming from F^. This cancellation 
is however not mathematically exact since it holds at the order 1/Ng, which indeed reduces to zero in the large sample 
limit. 

It is clear that the proof of such a tricky cancellation must rely on a very precise procedure in which the composite 
nature of the excitons is treated exactly. This is why the present work heavily uses the composite boson many-body 
theory for Frenkel excitons we have just constructed^'^. The main results of this theory which is very new, are recalled 
in Section II for the reader to possibly follow the rest of the argument. In Section III, we start with a simpler two-body 
problem, namely the calculation of the Hamiltonian mean value for N = 2 Frenkel excitons, while in Section IV, we 
turn to arbitrary N. This mean value is compared to a similar quantity obtained with excitons replaced by elementary 
bosons, (H) then only having a first order term in density. In Section V, we conclude. 




(1.6) 
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II. MAIN RESULTS OF THE COMPOSITE EXCITON MANY-BODY THEORY 



A. General results for composite excitons 



Many-body effects with composite excitons can be derived through a set of four commutatorsi^. The ones hnked 
to fermion exchanges are obtained from 



D ■ B^^ 

J-'im-) -'--'0 



(2.1) 



-7V(7V-l)i?r-2^Aao)st (2.2) 

n 

where definitions of the deviation- from-boson operator Dmi and the Pauh scattering A between two excitons 

follow from these equations taken for = 1. 

In the same way, many-body effects linked to fermion interactions are obtained froroi^ 



(2.3) 



H,Bl, 



= NB, 







EnB. 



(2.4) 



where definitions of the exciton creation potential V^^ and the interaction scattering between two excitons ^ ^"j"*;^ 
follow from these equations taken for = 1. 



B. Frenkel excitons 



In the case of Wannier excitons, the index i stands for (Q^, Vi) where is the exciton center of mass momentum 
and Vi a relative motion index. Due to the tight binding approximation which in some sense freezes electron on top 
of hole and makes Frenkel excitons constructed on electron and hole on the same atomic site, these excitons do not 
have relative motion index; so that they only are characterized by a center of mass momentum Q. Their creation 
operator Bq reads in terms of free pairs B^^ — al^bl^ on atomic site n as 



1 y ^-.Q.R„5t 



(2.5) 



where Ng is the number of these atomic sites and R„ their location. In the same way, free pairs read in terms of 
excitons as 



1 



as easy to check from the system periodicity which makes the following equations valid 

^eW-Q).R„=7V,5Q,Q 

n 

Q 

C. Fermion exchange 



(2.6) 

(2.7) 
(2.8) 



The fact that Frenkel excitons are only characterized by their center-of-mass momentum, leads to results linked to 
their many-body exchanges far simpler than the ones for Wannier excitons. In particular, the Pauli scattering of two 
Prenkel excitons is structureless: It only contains the expected momentum conservation between the "in" excitons 
(Qi, Q2) and the "out" excitons (Q'^, Q2) 



5, fQ'2 Q2\ 1 . 



(2.9) 



More generally, the Shiva diagram for carrier exchange between N Prenkel excitons shown in Fig. 1(b) leads to 

1 



N- 



(2.10) 



By using the Pauli scattering given in Eq. (2.9), we find that the commutators (2.1, 2) for fermion exchanges reduce 



to 



2N 



Q-Q' 



Bq,bI^ 



NBl--^ iS^o - Deo) - ^^^Bl--^Bl 



(2.11) 



(2.12) 



Prom them, it becomes easy to show that the factor Fn, defined in Eq. (1.3), which makes A'' Prenkel excitons differing 
from N elementary bosons, obeys the recursion relation 



"AT 



So that this factor is just given by 



Fn = 



Fn-1 



{Ns - N)\N^ 



(2.13) 



(2.14) 
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This compact expression of Fn as well as the rather simple commutators for fermion exchanges between N excitons 
given in Eq. (2.11, 12) are of great help to get exact results on Frenkel exciton many-body effects as shown below. 

D. Fermion interaction 

If we now turn to the interaction scattering induced by the Coulomb interaction which exists between electrons in 
the various atomic levels, we have been led to split it as 

(2.15) 

Actually, the last part ^neut of the interaction scattering ^ does not play any role in the many-body physics of 
Frenkel excitons because, when inserted in equations like Eq. (2.4), it gives 

E ^neut{$^^l)Bl,Bl,=0 (2.16) 

Q'iQ2 

Indeed, this neutrality scattering, which comes from the energy Eq we must pay to separate electron from hole, is 
equal to 

^neut (q? q') = 2£oA (^q? (2.17) 

while, due to Eqs. (2.5, 8) and the fact that B}^ = 0, we do have 

E ^ (S qO Bl,Bl, = ^ ^ = (2-18) 

QiQi " q 

The second part ^transf of the interaction scattering ^ comes from the indirect Coulomb processes between atomic 
levels, these processes insuring the excitation transfer over the whole system in the case of Frenkel excitons. Such 
a scattering, quite specific to Frenkel excitons, can actually be seen as a " transfer assisted exchange" , the coupling 
between the "in" excitons (Qi, Q2) and the "out" excitons (Qi, Q2) occurring through the Pauli scattering of these 
excitons 

^transf (g? = (Vq; + Vq^)A g^) (2.19) 

Its energy- like prefactor is related to the Q dependent part of the "out" exciton energy, namely 

Vq = E ^"'''■''^R (? o) (2.20) 



where Vr ^ is the Coulomb potential between atomic sites at R, one electron going from the atomic state to 
v[ while the other goes from V2 to , the relevant states for Frenkel excitons being the ground state v — Q and the 
first excited state v = 1. This Coulomb coupling Vr i^J, reads in terms of atomic wave functions as 



(:? = j dridr2^:,(ri)^:.(r2)^^:^-|^^^.,(r2)<^.,(ri) (2.21) 

Note that this transfer scattering £,transf which does not depend on the momentum transfer between the "in" 
and "out" excitons, but only on the momenta of "out" excitons, has similarity with the "photon assisted exchange 
scattering" we have recently identified between two polaritonsi^. 

The first part ^coui of the interaction scattering ^ in Eq. (2.15) comes from direct Coulomb processes between 
electrons, between holes and between electrons and holes. As for the Coulomb scatterings of two Wannier excitons, 
this scattering depends on the momentum transfer between the "in" excitons (Qi, Q2) and the "out" excitons (Q^, 
Q2) through 

4co«i (^q; qJ - ^Qi+Qi,Qi+Q2 (.^•^^i 

The prefactor Wq is similar to Vq but for all possible direct Coulomb processes 

Wq = E (i i) + (0 0) - Vr a 0) - Vr (? !)] (2.23) 

If we now construct exchange Coulomb scattering between two excitons along the standard procedure, namely 

r&i 3:)-Ea(S p:)^(p^ t) (2-24) 
P1P2 

we find from Eq. (2.9) and the fact that all scatterings conserve momentum, that 

Qi) ]^'^Qi+Qj!,Qi+Q2Xl^(Q?+Q qI) (2.25) 

Q 

It is then easy to see that, due to Eq. (2.8), the sum over Q reduces to zero for ^co«; as well as for S,transf since the 
i? = term is missing from the sums they contain. Consequently, 

(in (Q'2 Q2\ _ n _ f-in / Q2 Q2\ (n 

^transf ya'^ Qi J - - ^coul ^Q'^ J (Z.ZV) 

In the same way, if we construct the exchange scatterings for S,transf and ^coui between more than two excitons by 
adding a Coulomb line on the Shiva diagram shown on Fig. lb, we find from Eq. (2.10) that the Coulomb exchange 
scatterings shown in Fig. 2a in the case of four excitons reads as 



E ]v3'^'^'i+'^2+Q3+Q^P2+P3+Q3+Q4? (pi Qi) 



P1P2 



9 

= ]y5'^Q'i+Q2+Q3+Q4.Qi+Q2+Q3+Q4 (qi+Q Qi ) (2.27) 

* Q 

So that, this four-body exchange, as well as all high order exchanges, contains the same sum over Q which cancels. 
This shows that all Coulomb exchange scatterings between Prenkel excitons of the type of Fig. 2a reduce to zero. A 
way to understand the cancellation of these scatterings could be to say that Prenkel excitons are made of electron on 
top of hole; So that the "relative motion wave function" of these pairs reduces to zero. However this understanding 
is far too naive since Pauli scatterings for carrier exchange between two excitons do exist, as seen from the fact that 
Fn is not equal to 1 but exponentially small. 

Actually, there are two types of exchange Coulomb scatterings: we can either have the one of Fig. 2a in which 
all the "in" and "out" excitons are constructed on different pairs, but we can also have the one of Fig. 2b in which 
one exciton stays made with the same pair. This other Coulomb exchange scattering which, due to momentum 
conservation, reads as 

X] '^Q2+Q3+Q4.P2+Q3+Q4^ (q; Qi) 

P2 * 

_ 1 c , , , , f/Ql+Q2-Q'i Q2 \ /r,r,nx 

~ ]Y2"Qi+Q2+Q3+Qi.Qi+Q2+Q3+Q4? [q[ qJ [Z.Zb) 

does not contain the sum over Q appearing in Eqs. (2.25) or (2.27); so that it does not reduce to zero. 

Before going further, it can be of interest to note that Coulomb exchange processes like the one of Pig. 2c a priori 
exist for Wannier excitons. However, due to the symmetry between electron and hole in the case of Wannier excitons, 
the 2x2 direct Coulomb scattering of these excitons is such that 

?(r/)=0 (2.29) 

So that Coulomb exchange scattering like the ones of Pig. 2c reduces to zero for m = i. In contrast, since the wave 
functions for electrons in the atomic ground state and first excited state which enter the definition of ^coui through 
(4 1^1) defined in Eq. (2.21) are different, such an electron-hole symmetry does not exist for Prenkel excitons. 
This makes scattering like ^coui (q? q^) different from zero for Q'^ = Qi. 

At this stage, we can roughly say that exchange Coulomb scatterings like the ones of Pig. 2a are the relevant ones 
for Wannier excitons, as they are responsible for the r] expansion of (-ff) jv while this type of scatterings reduces to zero 
for Prenkel excitons. In contrast, diagrams like the ones of Pig. 2b are going to control Coulomb exchange processes 
for Prenkel excitons while they reduce to zero for Wannier excitons when one of the "in" excitons stay in the same 
state. 
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All this tends to show that fermion exchanges between Frenkel excitons are unimportant in some configurations 
but important in a few others. So that a careful calculation of these exchange processes through a secure many-body 
theory as the one we have developed, is highly desirable to possibly trust the obtained results. 

III. HAMILTONIAN MEAN VALUE FOR N = 2 EXCITONS 

As for many other problems we have already tackled, it is quite valuable to first perform calculations fov N = 2 
because these calculations are far simpler while most of the many-body physics already shows up. 
In order to calculate {H)j^ defined in Eq. (1.2) for N = 2, we first note that, due to Eqs. (1.3), (2.3) 

= 2(l-i-) (3.1) 

while Eq. (2.4) taken for N = 2, leads to 

HBl' \v) = 2EoBi' \v)+Y,^ °) B^Bl^ \v) (3.2) 

Q 

since Vq \v) = as easy to check by taking Eq. (2.4) for = 1 acting on |v).The Hamiltonian matrix element between 
two ground state excitons thus appears as 

{v\ Bl {H - 2Eo) Bl' \v}=Y,^ (q*^ i) {v\ B^B^bI^ \v} (3.3) 

Q 

To go further, we use the scalar product of two Frenkel exciton states which, according to Eqs. (2.1, 2, 9) reads as 

(^;| Bq,Bq,b1^^bI^^ \v) = [Sq, q,6q, q, - a gj)) + (Q^ ^ Q^) (3.4) 

When inserted into Eq. (3.3), the exchange Coulomb scattering generated by the A term of this scalar product reduces 
to zero due to Eqs. (2.24) and (2.26). We are thus left with 

{H), = 2Eo + ^^^^ (3.5) 

In view of this result, we are led to expect the Hamiltonian expectation value for A'^ excitons to have a N{N — l)/2 
prefactor in front of the interaction term ^ (q q) of {H)^ since the two excitons involved in this scattering are now 
taken among N; plus possibly higher order terms in N{N — 1){N — 2)..., due to fermion exchanges between more 
than two excitons. In addition, the term 1/A^s in the denominator of Eq. (3.5) could possibly be replaced by 

{N — 1)/Ns — ... (A — 1)(A — 2)/Ng + ... due again to fermion exchanges appearing in Fjv. 
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We are going to show that the amount of ./V-body exchanges in the Coulomb term exactly cancels the exchanges in 
Fn; so that the interaction term of {H)^ ends by reading as the one for {H)2 with just a N{N — l)/2 prefactor and 
nothing else. Let us now prove this unexpected cancellation which actually holds, within terms in l/Ns- 

IV. HAMILTONIAN EXPECTATION VALUE FOR N FRENKEL EXCITONS 

According to Eq. (1.3), the denominator of {H)^ in Eq. (1.2) is just NIF^. To calculate the numerator, we use 



Eq. (2.4) for 



tJV 



H, Bg 



. This readily gives 



+ 



N{N-1) 



QiQ; 



f /'Q2 o\ 
^ VQi 



NIFn 



V) 



(4.1) 



A. Algebraic derivation 



Since scattering conserves momentum, ^ differs from zero for Qi + Q2 = only. So that to get (-ff)jv) we 

must calculate 



Gn{Q) = {v\B^b1^b1^bI''-^\v) 



(4.2) 



As for Q = 0, this scalar product is just NlFj^, we in fact need to calculate Gjv(Q) for Q ^ only. Since Dqo \v) = 
as seen from Eq. (2.12) taken for A'' = 1 acting on vacuum, this Eq. (2.12) readily gives 



G*AQ ^ 0) = -^^^^^ (.1 B^'B^Bi^Bi^-^ \v) 



To go further, we use 



Bq,bI^ 



1 — -Dqq. This leads to split the RHS of the above equation as 



G^(Q^ 0) 



N{N - 1) 



{v\Br'Bl^-'\v) 



(4.3) 



-{v\B---^DcciBi''-'\v) 



+ {v\Br'Bl,BQBi''~'\v) 



(4.4) 



The first term is just {N - 2)!iijv-2- The second term, calculated from Eq. (2.11), gives [2(iV - 2)/iVs] {N - 
2)!F]v-2- In the last term, we use Eq. (2.12) for BqBq^~^ \v}. This allows us to rewrite this term as 
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— [{N — 2){N — 3)/Ns] Gjv-2(Q)- By collecting these three terms, we get a recursion relation between the Gjv's. 
It reads 



iV(Ar-l)(7V-2)(7V-3) 
+ Giv-2(Q) (4.5) 



Using Eq. (2.13) for Fjv, it is then easy to check that this recursion relation is fulfilled for 



Gn{Q^O) = -j;^N\Fn (4.6) 



So that we end with a surprisingly simple expression for G^iQ), namely 



When inserted into the sum of Eq. (4.1), we find that this sum reduces to 

o\ Ns6q 0-1 

Q 

Since the s\im of Coulomb scatterings over Q gives zero according to Eqs. (2.25, 26), we get the following compact 
form for the Hamiltonian mean value of N Frenkel excitons 

It shows that the interaction term is indeed the one for two Frenkel excitons with just a N{N — l)/2 prefactor which 
corresponds to the number of ways to choose these 2 excitons among A''. 

B. Shiva diagram derivation 

Although the above algebraic derivation is quite easy to follow, it may leave the reader unsatisfied as it does not 
really show why the amount of exchange processes in the Coulomb term is exactly the same as the one appearing in 
the normalization factor within terms in 1 /Ns . To grasp this reason, Shiva diagrams once more are quite valuable. 

Let us consider the Coulomb term appearing in {H)j^, namely 

^- = -T^ E i-\B^«Bi--^M^{to) (4.10) 
* Q1Q2 

In it, two excitons have a direct Coulomb interaction to become (Qi, Q2). They can then possibly mix with the 
other (N — 2) excitons through fermion exchanges to end as A'' excitons 0. This sum thus corresponds to to the 
diagram of Fig. 3a. 
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The simplest of these exchanges corresponds to leave {N — 1) excitons unaffected. So that these {N — 2) excitons 
appear as {v\ Bq~'^Bq^~'^ \v) = {N — 2)!F;v-2- The two other excitons (Qi, Q2) become (0, 0) with or without 
exchanging their fermions. This first possibility which corresponds to the diagram of Fig. 3b, thus leads to 



the N{N — 1) factor in the second bracket comes from the number of ways to choose the two excitons on the left 
among N. 

A second possibility is to mix one exciton with the excitons (Qi, Q2) while leaving the other {N — 3) excitons 
apart. This second possibility corresponds to the diagram of Fig. 3c. Wc have shown in Eq. (2.27) that Coulomb 
exchange scatterings like the one of Fig. 2a reduce to zero. By using the value of the exchange Shiva diagram of Fig. 
lb given in Eq. (2.10), we then find from the diagram of Fig. 3c that this second possibility produces a term 



in the Shiva expansion of Cn- The factor N{N — 1)(A^ — 2) comes from the number of ways to choose the three 
excitons of the left among N while the factor {N — 2) comes from the number of ways to choose the exciton on the 
right among {N — 2). This term has an overall minus sign because one fermion exchange is involved. An additional 
factor 2 comes from the number of diagrams which give the same nonzero contributions and which have the structure 

like the diagram of Fig. 2b. 

A third possibility is to mix two excitons with the excitons (Qi, Q2). The nonzero Coulomb exchange diagrams 
are shown in Fig. 4. They leads to 



[{N - 2)!F;v-2] [N{N - 1) (g S) - r (8 S))] 



(4.11) 



[{N - 3)!F;v-3] \n{N - 1){N - 2)-l$ (° °) {N - 2)1 2 



(4.12) 



+ [{N-4)\Fn_,] iV(iV-l)(iV-2)(iV-3)— ^(°°)(iV-2)(iV-3) 3 



(4.13) 



and so on... 



This Shiva diagram expansion of the Coulomb term Cat defined in Eq. (4.10) allows us to write it as 



Cn = Sn^ (0 0) 



(4.14) 



Sm = {N- 2)\Fn-2N{N -1)-{N- 3)\Fn-3N{N - 1){N - 2)-^{N - 2) 
+ {N- 4y.FN-4N{N - 1){N - 2){N - 3) (o 0^ (jv - 2){N - 3) - ... 



(4.15) 



The last step is to relate Sn to Fn- A simple way to do it is to consider the scalar product 



Pn{Qi, Q2) = (^;| B^Bl^B^^Bl''-^ \v} 



(4.16) 
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shown in Fig. 5a. Let us perform a Shiva diagram expansion of this scalar product according to the standard rules^^, 
namely we isolate {N — 2), {N — 3), ... excitons and write all the possible diagrams in which the remaining excitons 
are "not alone". This leads to the diagram of Fig. 5c for the second term. And so on...: By using the results for 
Shiva diagrams in the case of Frenkel excitons given in Eqs. (2.9, 10), we get for Qi = Q2 = 0, i.e., when PAr(Qi, Q2) 
reduces to NIFn 

NIFn = (iV - 2)IFn^2N{N - 1) (1 - 
+ {N - 3)!F^_3iV(iV - 1){N - 2) (^^ - Jj) (A^ - 2) + ... 

= (l-i-)5. (4.17) 

When inserted into Eqs. (4.14) and (4.1), this result allows one to readily recover the expression of the Hamiltonian 
mean value for N Frenkel excitons given in Eq. (4.9). 

It is clear that such a derivation needs to be knowledgeable with Shiva diagrams. It however has the great advantage 
to pick out in a transparent way from where the final result comes. In the present case, this is from Coulomb exchange 
processes like the one of Fig. 2. 



C. Physical understanding 

In order to identify the microscopic Coulomb processes which control the ground state energy of N Frenkel excitons, 
we first use the definition of the interaction scattering given in Eqs. (2.15, 16, 19). This leads to 

N(N - I) 

{H)^^NEo + ^^j^—^{Wo~2Vo) (4.18) 

When compared to Eq. (1.4), the microscopic value of the Coulomb contribution to the energy is thus found to read 
70 - Wo - 2Vo. 

To get a precise understanding of the various contributions to the Hamiltonian mean value, let us now write the 
above quantities in terms of the various elementary Coulomb scatterings between atomic sites. For that, we use Eqs. 
(2.16, 17, 18, 19, 20) and the expression of the Q exciton energy, namely Eq — Epair + Vq where the free pair energy 
Epair = Se + £h ~ £o with Eq = Vq (q J) — Vq (" J) is just the energy of one electron in the excited level and one hole 
in the ground state level, minus the neutrality energy Eq we would pay for not having electron and hole on the same 
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atomic site (see Reir^ ). This leads to 



{H)n - N 



R 



(4.19) 



since in the last sum, we can replace 2Vr, (° J) by Vr (" q) + ^ (o i)' *o ^'i- (2-21)- 

We see that the interaction term, which cancels for = 1, as expected, contains all possible Coulomb processes 
between the atomic levels and 1 of the different atomic sites, through the sum over R 7^ 0. Due to the orhogonality 
of the ly = (0, 1) atomic states, it is of importance to note that the indirect Coulomb processes Vr (" q) in which 
electrons change from atomic state 1 to atomic state are expected to be much smaller than the direct Coulomb 
processes Vr (J J) ; so that the interaction term of the Hamiltonian expectation value is essentially controlled by direct 
Coulomb processes Vr [^^^ . Let us however stress that, while these indirect scatterings are not that important in 
the interaction term of the Hamiltonian mean value, they play a major role in the Frenkel exciton physics because, 
although extremely small, they are the only processes allowing an excitation transfer between sites as necessary to 
produce the exciton. 

Since the Coulomb matrix elements Vr, [^^^ entering the Hamiltonian mean value are defined in terms of wave 
functions for electron in the ground and first excited atomic states, they strongly depend on the chemical properties 
of the material at hand. This has to be contrasted with Wannier excitons which are constructed on free electrons 
and free holes, the material of interest just appearing through effective masses induced by the lattice periodicity. For 
these Wannier excitons, the intraband Coulomb matrix element thus reduces to Ane'^/L^q'^ in 3D. As a result, the 
Hamiltonian mean value {H)j^ for these excitons only depends on one input parameter which is the exciton Bohr 
radius; this is why its i] expansion can be derived explicitly for all materials (see Eq. (1.6)). In contrast, the prefactor 
of the rj term for Frenkel excitons highly depends on the material at hand. 

If we now come back to Eq. (4.18), we can, in the large sample limit, replace Ng — 1 hy Ng. So that, for N large 
compared to 1, the Hamiltonian expectation value expands as given in Eq. (1.5): It only contains a linear term in 
density. Again, this has to be contrasted with {H)j^ for Wannier excitons in which the Coulomb term reads as a 
density expansion in rj. 

For completeness, let us note that the Hamiltonian mean value for N identical elementary bosons also has just 
one linear contribution in density. This is expectable since the Hamiltonian mean value contains the interaction at 
first order only by construction: when the particles are taken as elementary bosons the 2x2 interaction couples 
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two bosons among N, no fermion exchange exists to possibly mix these two particles with the rest of the system to 
produce higher order couplings. This is supported by the precise calculation of the Hamiltonian mean value 

, . {v\B^HBr\v} 
{v\BoBo \v) 



which for bosonized excitons such that 
effective bosonic Hamiltonian written as 



Bq',Bq 



<5q'Q) has a denominator which reduces to N\ while, for an 



H = ^EciBIBq + ^^V (^l ^l) B^q,B^q,Bq,Bci, (4.21) 



its numerator follows, since BqBq \v) = N6q qBq \v), from 



(.1 < (H - NEo) Bl- \v) = (3? I) 



{v\B^Bl,Bl,Blf'-'\v) (4.22) 
This leads to an Hamiltonian mean value given by 

(H)^ = NEo + ^^^^^V{ro) (4.23) 

By noting that the scattering V (o 8) must decrease as with the sample volume, as necessary for the extensivity 

of the system, we see that the correction to the energy of one bosonized exciton just increases linearly with the exciton 
density N/L'^. 



V. CONCLUSION 

In this paper, we calculate the ground state energy of N Frenkel excitons in the Born approximation through the 
Hamiltonian mean value in the state made of N ground state excitons. We show that the energy change induced 
by exciton interactions contains one linear term only in the exciton density, while higher order terms appear in the 
case of Wannier excitons. A naive reason could be that, since Frenkel excitons are constructed on highly localized 
atomic wave functions, their exchange Coulomb scatterings must reduce to zero, making Frenkel excitons appearing as 
elementary bosons. This however does not hold since the normalization factor of N Frenkel excitons is exponentially 
small compared to the one for N elementary bosons. The deep nonobvious reason for having one linear term only in 
the Hamiltonian mean value, is the fact that the same exponentially small prefactor also exists in the Coulomb term 
of the Hamiltonian matrix element, within corrections in 1/Ng where Ng is the number of atomic sites. So that, in 
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the end, all the exchange processes appearing in (H)^ finally disappear, as if Frenkel excitons were true elementary 
particles, within terms in l/Ng. Here again, Shiva diagrams enlighten the physical understanding of the fermion 
exchanges which lead to this rather unexpected cancellation. 
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Figure 1: (a) Shiva diagrams for carrier exchanges between two Prenkel excitons. This diagram represents the Pauli scattering 
A given in Eq. (2.9). (b) Shiva diagram for carrier exchanges between N Prenkel excitons, as given in Eq. (2.10). 
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Figure 2: (a) Coulomb exchange processes between 4 Frenkel excitons in which all "in" excitons and "out" excitons are made 
with different electron- hole pairs. This type of processes reduces to zero for Prenkel excitons (see Eq. (2.27)). (b) Coulomb 
exchange process in which one Prenkel exciton keeps its carrier pair (see Eq. (2.28)). (c) Same as (b) for Wannier excitons 
characterized by i = (Qiji^i). This type of Coulomb exchange process reduces to zero for Wannier excitons when the exciton 
which keeps its fermion stays in the same state {m = i). 
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Figure 3; (a) Shiva diagram representation of the Coulomb part of the Hamiltonian mean value, as given in Eq. (4.9). (b) First 
order term: (A'^ — 2) excitons stay unaffected in this process, (c) Second order term: one exciton among (A^ — 2) exchanges 
its fermion with the ones scattered by the 2x2 Coulomb interaction. 
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Figure 4; Non-zero Coulomb exchange process in which two excitons among {N — 2) exchange their fermions with the ones 
scattered by the 2x2 Coulomb interaction. 
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Figure 5: (a) Shiva diagrams representation of the scalar product of N exciton state in which 2 excitons on one side have 
different momenta Qi and Q2. (b) First order term: {N — 2) excitons stay unaffected. The excitons (Qi, Q2) can either 
keep their fermion or have an exchange, (c) Second order term: {N — 3) excitons stay unaffected. 



